Kondo Force in Shuttling Devices: Dynamical Probe for a Kondo Cloud 
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We consider electromechanical properties of a single-electronic device consisting of movable quan- 
tum dot attached to a vibrating cantilever, forming a tunnel contact with a non-movable source elec- 
trode. We show that the resonance Kondo tunneling of electrons amplify exponentially the strength 
of nanoelectromechanical (NEM) coupling in such device and makes the latter to be insensitive to 
mesoscopic fluctuations of electronic levels in a nano-dot. It is also shown that the study of Kondo- 
NEM phenomenon provides an additional (as compared with a standard conductance measurements 
in a non-mechanical device) information on retardation effects in formation of many-particle cloud 
accompanied the Kondo tunneling. A possibility for superhigh tunability of mechanical dissipation 
as well as supersensitive detection of mechanical displacement is demonstrated. 



PACS numbers: 73.23.-b, 72.10.Fk, 73.23.Hk, 85.85.+j 

Recent progress in fabrication of nanoelectromechan- 
ical systems (NEMS) based on suspended carbon nan- 
otubes [I] as well as on suspended Si [2 and SiN [3] 
nanowires vibrating at radio frequencies (RF) resulted 
in rapidly growing amount of theoretical works [J- [El 
addressing the issues of interplay between spin/charge 
transport and nano- mechanics [9 . The observation of 
Coulomb Blockade [2 in NEMS opened a possibility to 
consider the influence of strongly correlated and reso- 
nance effects on a behaviour of nano-oscillators. 

Usually, NEM regime implies strong coupling between 
the electronic and mechanical degrees of freedom. The 
coupling is provided by two main mechanisms. Motion 
of the movable dot (shuttle) between two metallic banks 
results in the time dependent tunneling amplitudes. On 
the other hand, the electron charge transport between the 
banks in presence of magnetic field results in appearance 
of a Lorentz/Laplace force acting on the shuttle, which 
should also be taken into account [4j [5] . The aggregate 
dynamics of a shuttle is that of a periodic oscillator with 
decrement or increment and an electron tunneling (co- 
tunneling) parametrically dependent on this slow classi- 
cal motion. In some sense the problem may be treated 
as a tunneling through anharmonic vibronic system. In 
many cases, e.g. in shuttling devices including bending 
carbon nanotubes [5 , the vibronic language completely 
describes the physical situation. 

NEM coupling like other nanometer length scale phe- 
nomena is strongly affected by mesoscopic fluctuations. 
Spatial quantization of electronic motion in a quantum 
dot makes electro-mechanical transduction to be sample 
sensitive phenomenon [2 . An exception to this rule is 
electromechanical coupling due to the many-body Kondo 
tunneling. Indeed in this case the charge transfer is 
controlled by the singularity of the tunneling density of 
states at the energy pinned to the Fermi level of the injec- 



tor and thus protected against mesoscopic fluctuations. 
This fact in combination with another generic feature of 
Kondo phenomenon - its super-sensitivity to a strength of 
the tunneling coupling (and therefore its super-sensitivity 
to the mechanical displacement of a quantum dot) - 
makes Kondo NEM coupling to be a phenomenon promis- 
ing for practical applications. 

An example of such device is schematically shown in 
Fig. [I] A nanoisland is mounted on the metallic can- 
tilever, which may vibrate under an external force. The 
contact between the source and drain electrodes is a com- 
bination of time-dependent tunneling bridge between the 
source S and the island and a metallic bridge formed by 
vibrating cantilever connecting the island with the drain 
D (see Ref. [2] for experimental realization). 

We consider the configuration where a cantilever is dis- 
places in y direction, u = (0, 0) in a magnetic field B = 
(0, 0, B). In this case the Laplace force F acts on the can- 




FIG. 1. (Color online) Shuttle with a cantilever. The nanois- 
land (green) is separated from the wire (drain) by a tunnel 
barrier (light blue) of constant width. The width of the bar- 
rier to the source changes during the cycle. Shuttle oscilla- 
tions are stimulated by initial conditions, e.g. bending the 
cantilever by the gate voltage at time t = 0. 
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tilever in the same direction y, Fl = L ■ I x B = (0, F, 0). 
Here L is the length of the cantilever. Besides, small elec- 
tromotive force F em f = (/, 0, 0) acts on the electrons in 
the cantilever. In the limit of strong Coulomb blockade 
in the nanoisland, the Kondo screening accompanies the 
tunneling 'source-island-cantilever', and a unique possi- 
bility arises to study the contribution of purely quantum 
many-particle Kondo effect on the classical oscillation of 
a shuttle (cantilever + island). The study of a "Kondo 
force" in shuttling is the main subject of this paper. 

We study two coupled subsystems: a tunneling con- 
tact 'source-moving island-moving cantilever' treated as 
a purely quantum system in a framework of the Anderson 

- Kondo model, and a macroscopic wire with attached 
island oscillating under external constraining force. We 
work in the Kondo limit, where the nanoisland is repre- 
sented by its spin 5, so that internal degrees of freedom 
are the spin-flip processes. The source-drain transport 
is a combination of quantum tunneling 'source - island 

- cantilever' and ohmic transport 'island-drain'. In this 
"time-dependent Schrieffer - Wolff' limit (see below) the 
Hamiltonian of quantum subsystem is 
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Here the indices /,r stand for the electronic states in 
the source and cantilever, respectively, = e k — [i are 
the excitation energies of lead electrons, N a — ^2 ka n aka 
are the corresponding electron density operators, s aa , — 

operators for the electrons in the leads and in the nanois- 
land, respectively, r and a are the vectors of Pauli ma- 
trices acting on the states in the leads and dot. At small 
bias voltage Vbias the source and the cantilever are sup- 
posed to be in the adiabatically stationary state of ther- 
mal equilibrium. The parameters are J a0i > = 4v^v a , / E c ^ 
where v a is the tunneling amplitude between the nanois- 
land and the metallic lead a, E c is the Coulomb blockade 
energy. The exchange couplings J u and J lr are time- 
dependent due to the dependence of the tunneling am- 
plitude between the source and the moving nanoisland on 
the island position v l = ^ [£?(£)]. The time dependence 
of this amplitude is a set of pulses corresponding to elec- 
tron injection from the metallic reservoir to the shuttle 
periodically approaching the bank S. [10] We confine our 
treatment with the simplest case of S = 1/2 (odd occu- 
pation of a nanoisland in the neutral state) and the single 
channel tunneling between the nanoisland and the leads. 

The oscillations of cantilever with attached nanoisland 
are determined by the classical Newton equations 



— u + UJnU = —F. 
Qo m 



(2) 



where ujq = ^k/m is the oscillator frequency of free can- 
tilever, Qo is a quality factor of NEM device. 

Our aim is to study the spin and charge transport by 
means of a shuttle oscillating in accordance with Eq. ([2| 
in presence of many-particle Kondo screening described 
by the Hamiltonian ([!]). The coupling between the clas- 
sical and quantum subsystem is realized via the param- 
eters Ju(u), Jir(u), F(u), where the time dependence 
u(t) should be calculated self-consistently. Meanwhile, 
J rr does not depend on displacement u (see Fig. 1). 

The cotunneling Hamiltonian may be rationalized by 
means of the Glazman-Raikh rotation, which in our sit- 
uation is time dependent: 
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with tan$£ = \v r /vi(t)\. Use of this transformation 
for diagonalization of the Schrodinger operator C = 
—ihd/dt + H(t) results in generation of additional term 
Hb proportional to —ihU^ 1 d t U in the transformed 
Hamiltonian (see, e.g., [14]). 
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The term Hb may be treated as an additional gauge 
potential in the lead Hamiltonian describing a Berry-like 
phase [14] generated by shuttle motion [13]. Only the 
even partial wave ipi survives in the cotunneling term 
H ex with the time-dependent effective indirect exchange 
coupling J(t) = Ju(t) + J rr (see, e.g., [IS]). This time 
dependence may be parametrized [9] in assumption that 
the source-island tunneling amplitude is an exponential 
function of a distance y between the source and the 
moving nanoisland, while the tunneling nanoisland- 
cantilever is constant: v r = vq, vi = v\ exp[?/(£)/A]. 
The spatial coordinates are counted off the equilibrium 
position of cantilever, so that v\ is exponentially small, 
v\jv§ ~ exp(— yo/X). Here A is the confinement radius 
(tunnel length) of the electron wave function within the 
island, yo is the distance between the source and the 
island at equilibrium. 

We suppose that the shuttling mode is slow enough and 
the electron transport is adiabatic, i.e., the exchange cou- 
plings Ju , Jir depend parametrically on time via the dis- 
placement coordinates u(t). Then to find the tunneling 
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current one may trace the time dependence of local oc- 
cupations of the left and right banks (source and nanois- 
land) near the point of tunneling contact. The current 
operator is 



e d 
2dt 



(Nr-Nt) 



(6) 



where in the case of immovable nanoisland only the even 
mode 1 contributes to the current. In our case both 
modes 1, 2 are involved in the tunneling transport due 
to the term Hb in Eq. Q. After the Glazman-Raikh 
rotation the current operator transforms into 



where 



a- d A d A 
X= dt Qt +dt qt > 
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Qt = ~\ Sin 2^ (flka^ka 



Here the operator Qt controls the time-dependent elec- 
tron occupation in the source lead, and the operator 
q t is responsible for all tunneling and cotunneling pro- 
cesses including admixture of odd components ^ 2/ccr to 
the tunneling charge transport induced by the gauge field 

H B .[MUM 

The time-dependent Glazman-Raikh angle defined by 
([9| results in adiabatic time dependence of the Breit- 
Wigner factor 
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Using the Friedel - Langreth sum rule [T7] . one may 
write 



Ni - N 2 



(10) 



where St = S^+S^ is a total time dependent Friedel phase. 
At the unitary limit 5^ = ±7r/2. 

We are interested in the Kondo effect contribution to 
the tunneling current. This contribution is characterized 
by the spin dependent scattering phase shift S cr (e) in the 
source lead, which approaches the unitarity limit ir/2 at 
T— ^0 and e^e F . In the adiabatic limit Huj <C k B T^ in 
under conditions (k B T, gfi B B, \eV hias \) < k B T™ in the 
phenomenological Fermi liquid Hamiltonian i^Noz may 
be used [18 (here k B is Boltzmann constant, \i B is 
Bohr magneton and g is Lande factor). In this Hamil- 
tonian both scattering and interaction are scaled by the 
time-dependent Kondo temperature Tx{t) taking mini- 
mal value 7™ n at maximal distance from the source. 

In order to get full tunnel current in adiabatic ap- 
proximation we (i) calculate a linear response with re- 
spect to both bias eVbias <^ k B Tx and hd$ t /dt < 



huiofimax <C k B Tx, (h) take into account cancellations 
arising due to emergent SU (2) symmetry associated with 
channels [T3].[20]. (hi) perform averaging with the adia- 
batic Hamiltonian (|4][5| at zero bias and zero temper- 
ature. The finite temperature and bias effects are ac- 
counted by Nozieres method. [18], [19] As a result, the 
tunnel current l t = Io(t) + I[ n t(i) consists of two parts: 
the Friedel phase contribution 
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and G^=i 2 are bold retarded (A = R), advanced (A = A) 
and Keldysh (A = K) Green's function [2T] , 

Let us first rewrite the Friedel part of the tunneling 



current ( 11 ) via the parameters characterizing the Kondo 



tunneling in the low-temperature strong coupling limit at 
T < T™ in , where 



k B T K (t) = D exp 
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Dq is the ultraviolet cut-off for the Kondo problem with 
a scale of the band energy in the source, T a = 7rpo\v a \ 2 , 
po is the density of electronic states at the Fermi level 
£p. In our adiabatic regime Tk parametrically depends 
on time, following the time dependence of Ti(t). The 
Hamiltonian i^Noz [IE], [19] establishes the relations be- 
tween 8 (j, B and Tk near the unitary limit, such as 
S t = 2\eV hias \/(k B T K (t)) < 1. The magnetic field en- 
ters only into relative Friedel phase 5^- — 5± = tt — 
2gp B B/(k B TK(t)). We neglect the influence of magnetic 
field on Tk, since we work in the limit gp B B <C 7™ n . 
Alternatively, a non-uniform magnetic field negligible at 
the dot and gradually increasing along the cantilever 
could be assumed in the model. 

In the adiabatic limit the Friedel phase S t and Glaz- 
man - Raikh angle $ t are not independent, but connected 
through (|9) and (13) 

1 dS t ttE c . dd t nA s 
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with To = 7rpo\vo\ 2 <C E c . Thus, the Friedel contribution 
to tunnel current can be expressed in terms of shuttle 
velocities as follows: 
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FIG. 2. (Color online) Time dependence of the cur- 
rent Io for different values of asymmetry parameter u = 
yo/X. For all three curves shuttle oscillates with ampli- 
tude ym ax = A, hjujo/(k B T% in ) = 10" 

gii B B/{k B Tr n ) = 0.1 with T*™ 
ln[Qo/Q] as a function of u, Qo 



" 3 , \eV hias \/(k B T^ n ) = 
2K, X/L = 10" 4 . Insert: 
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Here the time dependence of the tunnel current is prede- 
termined by the time dependence of tunnel integrals for 
the nanoisland moving in y direction, i.e. by the function 
y(t) and its derivative y. Moreover, one can see that even 
in the case of possible instability large amplitude oscilla- 
tions are exponentially suppressed. Typical behaviour of 
Io(t) is shown in Fig (2] The non-sinusoidal form of cur- 
rent is associated with time dependence of both tunnel 
width and Kondo temperature. 

The second term L m t(t) given by Eq. (12) leads to 



"ohmic" contribution to the current with unitary con- 
ductance Go = e 2 /h (see discussion below): 

lint W = G V hias sin 2 20 t Yl sin2 5 ° ( 16 ) 
The force in the r.h.s. of the Newton equation Q is 
a sum of the driving force F , the Laplace force Fl and 
electromotive (emf) force F em fi 



F(y,t) = F (t) +l t • B • L + F emf 



(17) 



The emf force can be estimated as F em f~y(B • L) 2 Go 
[22] . Due to sequential geometry of electric circuit, the 
current X t = Io(t) + Imt(t) is the tunneling current de- 
fined by ( pi)] ) and (16) . In the limit of small bias volt- 
age |eVbias| kBT^ 171 , electrons in the source and the 
cantilever are supposed to be in adiabatically station- 
ary state of thermal equilibrium. Then the parametriza- 



tion (15) is valid and with accuracy to small param- 
eters 0([(eV bias /(k B T™ n )] 2 , [(g^ B B)/(k B T^ in )] 2 ) the 
Lorentz force may be written as 



F L = F ad (y(t)) - y 



dF n , 



HnEr 



(18) 



where F ad = 2B • L • G • V h [ as cosh" 



-2 [y(t)-y ] 



and T\ 



(o) 



K 



is a Kondo temperature at equilibrium position. Small 



correction to the adiabatic Lorentz force in the (18) may 



be considered as a first term in the expansion over a 
small non adiabatic parameter ujqt <C 1, where r is the 



retardation time associated with inertia of the Kondo 
cloud. Using such interpretation one gets r: 



h7rE r 
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where Q(B) and uj(B) are the quality factor and oscil- 
lator's frequency at finite magnetic field B respectively. 
Equation Eq{l9] allows one to obtain information about 
dynamics of the Kondo clouds from the analysis of the 
experimental investigation of the mechanical vibrations. 
The retardation time associated with dynamics of Kondo 
cloud is parametrically large compared with the time of 
formation of the Kondo cloud tk = Ti/Tk and can be 
measured owing to small deviation from adiabaticity. [23] 
Also we would like to emphasize a supersensitivity of the 
quality factor to the change of the equilibrium position 
of the cantilever characterizing by the parameter yo . The 
plot ln[Qo/Q] is presented in insert of Figj2] From this 
plot one can see that both suppression Q > Qo and en- 
hancement Q < Qo of the dissipation of nanomechanical 
vibrations (depending on the direction of the magnetic 
field and the equilibrium position of the cantilever ) can 
be stimulated by Kondo tunneling. The latter demon- 
strate potentialities for the Kondo induced electrome- 
chanical instability which will be a subject for separate 
analysis. 

Equations ( 15|18 ) and ( 19 ) represent the central results 
of the Letter. On the one hand, we have shown that the 
electric current associated with the Kondo effect results 
in magnetic field dependent Q-factor allowing to fine- 
tune the nano- mechanical resonator. On the other hand, 
the non-ohmic part of the current provides an informa- 
tion about retardation effects related to the motion of 
the Kondo cloud. Thus, the measurement of the Kondo 
forces in Single Electron Transistor give a complementary 
to conductance measurements information. 

In conclusion we have shown that the Kondo phe- 
nomenon in single electron tunneling gives a very 
promising and efficient mechanism for electromechani- 
cal transduction on a nanometer length scale. Measur- 
ing of nanomechanical response on Kondo-transport in 
nanomechanical single-electronic device enables one to 
study kinetics of formation of Kondo-screening and of- 
fers a new approach for studying nonequilibrium Kondo 
phenomena. Kondo effect provides a possibility for super 
high tunability of the mechanical dissipation as well as 
super sensitive detection of mechanical displacement. 
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